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Monte Carlo (MC) simulations are used to study the dynamics of polymer translocation through a nanopore 
in the limit where the translocation rate is sufficiently slow that the polymer maintains a state of conforma- 
tional quasi-equilibrium. The system is modeled as a flexible hard-sphere chain that translocates through a 
cylindrical hole in a hard flat wall. In some calculations, the nanopore is connected at one end to a spherical 
cavity. Translocation times are measured directly using MC dynamics simulations. For sufficiently narrow 
pores, translocation is sufficiently slow that the mean translocation time scales with polymer length N ac- 
cording to (r) (X {N — iVp)^, where A^p is the average number of monomers in the nanopore; this scaling 
is an indication of a quasi-static regime in which polymer-nanopore friction dominates. We use a multiple- 
histogram method to calculate the variation of the free energy with Q, a coordinate used to quantify the 
degree of translocation. The free energy functions are used with the Fokker-Planck formalism to calculate 
translocation time distributions in the quasi-static regime. These calculations also require a friction coeffi- 
cient, characterized by a quantity N^g, the effective number of monomers whose dynamics are affected by the 
confinement of the nanopore. This was determined by fixing the mean of the theoretical distribution to that 
of the distribution obtained from MC dynamics simulations. The theoretical distributions are in excellent 
quantitative agreement with the distributions obtained directly by the MC dynamics simulations for physi- 
cally meaningful values of A'cff. The free energy functions for narrow-pore systems exhibit oscillations with 
an amplitude that is sensitive to the nanopore length. Generally, larger oscillation amplitudes correspond to 
longer translocation times. 



I. INTRODUCTION 



The translocation of a polymer through a narrow pas- 
sage from one space to another is a fundamental process 
with numerous applicationsji It is a key part of various 
biological phenomena, including viral genome transfer, 
RNA transport through the nuclear pore complex, pro- 
tein transport into mitochondria, and genome packing 
in bacteriophages)^ The pioneering work of Kasianow- 
icz et al. first demonstrated the feasibility of monitor- 
ing translocation of single molecules by measurement 
of the ionic current passing through nanometre-sized 
pores4i^ In this method, the potential difference that 
drives the ionic current also drives a charged macro- 
molecule through a nanopore, and the translocation re- 
sults in a measurable reduction in the current. Substan- 
tial effort has been devoted to further develop this exper- 
imental technique using both biological®" ^■^ and synthetic 
nanopores ji^"— This work has largely been motivated by 
the desire to develop a fast and inexpensive method of 
sequencing DNA.— "— The experiments have revealed a 
complex dependence of translocation behaviour on a va- 
riety of factors, including electric field strengthjii pH,'^^ 
and salt concentration. ^-2, An understanding of this grow- 
ing body of experimental results is currently incompleteji 

Numerous theoretical and computer simulation stud- 
ies of polymer translocation have been reported in re- 
cent years. Much of this work has been summarized 
and evaluated in a recent review. One notable the- 
oretical approach, pioneered by Sung and Park'^^ and 
Muthukuniar'^® employs the Fokker-Planck (FP) formal- 
ism. Here, the dynamics are governed by the free en- 



ergy, F(s), which is a function of a single coordinate, 
s, typically defined to be the number of monomers on 
one side the nanopore. The time-dependent probability 
distribution for this coordinate is obtained by solving a 
1-D Smoluchowski equation, for which the free energy 
function is required input. This solution can be used to 
calculate translocation time distributions and determine 
scaling relations for mean translocation times. The pre- 
dicted scaling relations are quantitatively consistent with 
results from experiments of DNA translocation through 
a-hemolysini— The FP formalism has been used to 
study translocation to various levels of approximation 
and for different confinement geometriesJ^i^"— A de- 
tailed review of the FP formalism and its application to 
polymer translocation is provided in the recent book by 
Muthukumar.^ 

The veracity of the free energy functions clearly de- 
termines the accuracy of the predictions for translo- 
cation dynamics in the FP approach. Most FP the- 
oretical studies employ analytical approximations of 
i^(s)iii2^i^i22i"^i^"^ In other studies, the free energies 
have been calculated numerically for specific molecular 
models that can also be employed in dynamics simu- 
lations. For example, exact enumeration of polymer 
configurations has been used to calculate a free energy 
landscape in order to elucidate the effects of polymer- 
nanopore attractions on translocation times for a lattice 
model system4ii^ Off-lattice models have been stud- 
ied using specialized Monte Carlo simulation methods. 
One recent example is the incremental gauge cell (ICC) 
method,'''^ which was used to calculate free energy pro- 
files for a flexible Lennard- Jones chain translocating out 
of a spherically conflned space with adsorbing wallsj— 
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Another efRcient MC method appropriate to calculate 
translocation free energies is the self-consistent multiple 
histogram (SCMH) technique.— Recently, we used the 
SCMH method to calculate free energy functions for flex- 
ible hard-sphere chains translocating through a cylindri- 
cal pore in a hard barrier 4^*^ Generally, the free energy 
curves follow trends that conform well to the predictions 
of a modified version of the standard analytical model 
described in Ref.[l|. However, one unexpected feature ob- 
served was an oscillation in the free energy with an ampli- 
tude that is very sensitive to the nanopore dimensions.^- 
While its relevance to the experimental context is un- 
clear, knowledge of such effects may be helpful for the 
interpretation of results from dynamics simulations and 
could be relevant to elucidating the sensitivity of the 
translocation behaviour on seemingly minor differences 
in the model that have been observed in some simulation 
studies^^>^ 

The validity of the FP approach to polymer transloca- 
tion is founded on the assumption that the portions of 
the polymer outside the pore remain in conformational 
equilibrium during the entire translocation process. This 
is expected to be the case in the limit of sufficiently 
strong pore friction for polymers of finite length. In 
this case the translocation time scales as N'^ for unbiased 
translocatio n^^'^^ and as N/ f for forced translocation)^^ 
where / is the driving force magnitude. However, it 
is well known that this assumption must break down 
for sufficiently long polymers42ii^ Computer simulation 
studies have generally yielded scaling exponents incon- 
sistent with the values for the quasi-static regime with 
strong pore friction42r— i^i^""— The values of the ex- 
ponents vary widely and appear to be sensitive to the 
details of the model, simulation method, and the length 
of the polymer. In addition, nonequilibrium conforma- 
tional behaviour has been observed directly] ^^'^"^i'^^'''^' '— 
Various theoretical models have been developed to ac- 
count for non-equilibrium effects and generally predict 
scaling exponent values greater than those for the quasi- 
equilibrium case, both for nondriven^^— and driveit^^— 
translocation. Clearly, the FP formalism is an inappro- 
priate approach for the conditions considered in most 
of these simulation studies. However, simulation models 
can easily be adjusted to either reduce the relaxation time 
of the chain or increase the nanopore friction, either of 
which should eventually lead to quasi-static conditions. 
This is illustrated by a recent simulation study of un- 
forced translocation that followed the former approach 
and measured a scaling exponent value consistent with 
FP predictions for the quasi-static regime at sufficiently 
low viscosity where the relaxation time of the chain is 
very short 

The purpose of this study is to examine the accuracy 
with which the FP formalism describes the translocation 
dynamics in the quasi-static regime. The system is mod- 
elled as freely-jointed hard-sphere chain that translocates 
through a cylindrical nanopore in a hard barrier. We 
employ free energy functions that were calculated using 



the SCMH method in our recent study. The calculated 
translocation time distributions are compared to those 
measured directly from MC dynamics calculations for the 
same model. Calculations were carried out for systems 
in the quasi-static regime with friction dominated by the 
nanopore. The regime boundary is determined by mea- 
surement of the scaling exponent in MC dynamics simula- 
tions. The friction coefficient used in the FP calculations 
is characterized by a parameter Ncs, the effective number 
of monomers whose dynamics are affected by nanopore 
confinement. We focus mainly on a system with a pla- 
nar barrier and semi-infinite cis and trans spaces. We 
also present some results for the cases where the pore is 
connected to a spherical cavity at one end; comparable 
systems have been the focus of a number of other theo- 
retical and computer simulation studiesj ^^'^^i'^°'^^i'^^'^^i^^ 
We investigate the effects on the translocation dynamics 
of varying the nanopore dimensions, the initial position of 
the polymer, and polymer length. Generally, the theoret- 
ical distributions are in excellent quantitative agreement 
with those measured from MC dynamics simulations for 
physically meaningful values of A^cff- 

The remainder of this article is organized as follows. 
In Section [Hi the model is described. In Section IIIIl 
we present a brief overview of the relevant results of the 
FP methodology and describe how the friction coefficient 
employed in the theory is calculated. In Section HVl we 
describe the MC algorithm used in the dynamics simula- 
tions and the details of the FP calculations. In SectionlVl 
the simulation results are presented, and the significance 
of the results is discussed in Section IVTl 



II. MODEL 

The polymer is modeled as a flexible chain of N hard 
spheres, each with a diameter of a. The pair potential 
for nonbonded monomers is thus 



Unhir) 



oo, r < a 
0, r > cr, 



(1) 



where r is the distance between the centres of the 
monomers. Pairs of bonded monomers interact with a 
potential 



Ub(f) = 



0, 



0.9(7 < r < 1.1(7 



oo, otherwise. 



(2) 



Thus, the bond length can fluctuate slightly about its 
average value, which is equal to the monomer diameter. 

The polymer undergoes translocation through a nar- 
row pore from one region to another. The nanopore 
is modeled as a cylindrical hole of length L and radius 
i? in a hard barrier. The regions on the two sides of 
the nanopore are labeled cis and trans. Two different 
cases are considered in this study. In the first case, the 
cis and trans sides are semi-inflnite spaces, and the bar- 
rier that separates them is bounded by two flat surfaces. 



3 



Monomers cannot overlap with the barrier, and thus the 
minimum distance between the centre of a monomer and 
the nearest point on either the barrier wall surface or the 
nanopore surface is (t/2. This confinement geometry is 
illustrated in Fig. [Tfa). In the second case, transloca- 
tion occurs between a spherical cavity of radius Rc on 
the cis side to a semi-infinite space on the trans side, as 
shown in Fig. [Tib). Most of the calculations in this study 
employed the first variant of the model. 




FIG. 1. Illustration of the two confinement geometries em- 
ployed in this study. Most of the calculations in this study 
used the system shown in (a). 

The degree to which the polymer has translocated 
across the nanopore is quantified using a translocation 
coordinate, Q, which is defined: 

1 ^ 

Q^nT.Q^^ (3) 

1=1 



where the translocation coordinate for the ith monomer 
is given by 





'0,zj < 




g. - < 


zjL,0< z, <L 


(4) 




l,z, > L 


1 



where Zi is the z-coordinate for the ith monomer. Thus, 
(5 = when all monomers are on the cis side of the 
pore, Q — 1 when the polymer is completely on the 
trans side, and < Q < 1 if there are any monomers 
inside the nanopore. Clearly, Q is a continuous quantity, 
which makes it a more convenient choice for a translo- 
cation coordinate here than the quantity s (the number 
of monomers on the trans side of the pore) used in most 
translocation studies. 
III. THEORETICAL BACKGROUND 

In the quasi-static limit, the time-dependent translo- 
cation probability distribution, V{Q,t), is governed by 
the Fokker-Planck equation. When the pore friction is 
dominant contribution to the total polymer friction the 
form of the equation iai 



where fee is Boltzmann's constant, T is temperature, 
F{Q) is the equilibrium translocation free energy func- 
tion and D is an effective diffusion coefficient. Consider 
a domain bounded by Q = Qa and Q = Qb- The first 
passage time, r, for a polymer located at Qo G [Qa,Qh] 
at time t = to reach either boundary has a probability 
distribution given byi 



9iT;Qo) = -^ J^' dQpiQ,T;Qo,0), (6) 



where p{Q,t;Q',t') is the conditional probability of the 
polymer reaching a value Q at time t given that it had a 
value of Q' at an earlier time t'. For a Markov process, 
p{Q, t; Q', t') also satisfies Eq. ([5]). In the case of absorb- 
ing boundary conditions the mean first passage time is 
given byji 



/ rQb dy \ rQo dy' TT(.,f\ _ ( fQo dy \ rQb dy' TT(.,f\ 
\JQa ■>P{y))jQ. xp{y')^yy I \JQ^ ^{y) J JQo V(y') ^ 



= jQ^y ' 



I 

where Finally, the probability a polymer reaches Q — Qh with- 

V'(Q) = exp (-[F(Q) - FiQ,)]/kBT) . (8) ^^^'^^^'^S S^^^^ ^^'"^ 



and 



HiQ) = r dQ'i,{Q') (9) ^,(go) = l^g^ (10) 
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where 

Note that this probabihty depends only on the free en- 
ergy function and not on the diffusion coefficient. 

In order to calculate g{T;Qo) and (r), the free energy 
function, F{Q), and the diffusion coefficient D must be 
suppHed. In this study, F{Q) is calculated using the 
SCMH method, as described in Ref. 13 The diffusion 
coefficient D is closely related to the diffusion coefficient 
of monomers confined by the nanopore. In this study, 
we employ a MC method to study the translocation dy- 
namics of the polymer; the details of the method will be 
described below in Section IIV Bl The diffusion rate for 
monomers inside the nanopore is determined by nature 
of the MC trial moves for the individual monomers. Un- 
fortunately, quantitatively accurate values of D cannot 
easily be calculated analytically from knowledge of the 
details of the model and method alone. Consequently, 
the diffusion coefficient is essentially treated as a free 
parameter. However, it is possible to determine the scal- 
ing of D with polymer length and (approximately) with 
nanopore length in the limit of strong pore friction. 

Consider a polymer that moves inside an infinitely long 
cylindrical tube aligned along the z axis. The polymer 
diffuses along z such that its mean square displacement 
satisfies 

(Az2) ^ 2D,t = (12) 

where Dz — k^T/jz is the diffusion coefficient and 72 is 
the friction coefficient associated with the polymer inside 
the tube. In this study, we consider the model polymer 
described in Section HIl that moves via random displace- 
ments of randomly chosen monomers, one at a time. In 
this case, the friction coefficient in Eq. ([T2|) scales as 

7. =7o(A^-l), (13) 

where N is the number of monomers. In addition, 70 is a 
constant that depends on the radius of the tube and the 
probability distribution of trial monomer displacements 
employed in the calculations. This scaling relation is con- 
firmed through results from simulations. Now consider a 
polymer undergoing translocation through a cylindrical 
nanopore of finite length, L. For a sufficiently narrow 
pore, the pore friction is expected to be the dominant 
contribution to the overall friction of the polymer. In 
this case, the friction coefficient governing the motion 
along z for the portion of the polymer inside the pore is 
expected to obey a relation similar Eq. (|13l) . 

7. =7o(iVcff-l), (14) 

where N^g is the effective number of monomers whose 
dynamics are strongly effected by the presence of the 
nanopore. To estimate A^off, first note that the range 
over which monomers are subject to some degree of con- 
finement is L + a. For a narrow nanopore, the poly- 
mer is approximately linear inside the pore and the 



number of bonds spanning the confinement region is 
ribond = (L + cr) / a] thus, the number of monomers span- 
ning this region can be estimated as A^oft ~ "bond + 1, 
or 

N,s~L/a + 2. (15) 

In practice, N^g is a free parameter, albeit one that 
should come close to satisfying the approximation in 
Eq. (USD. 

In the absence of a free energy gradient, and for suf- 
ficiently strong pore friction, the mean square displace- 
ment for monomers inside the nanopore is expected to 
obey Eq. p^ . From the definition of Q in Eqs. ^ 
and (|4]) it follows that a pore monomer displacement of 
Az corresponds to a translocation coordinate change of 
AQ = Az/{N(j), on average. Consequently, (AQ^) = 
2Dt, where 

D^^-^, (16) 
7 

and where 7 = N'^jz- D is the diffusion coefficient that 
appears in Eq. ([5|), and 7 is the corresponding friction 
coefficient. It follows that 

7 = iV27o(JVeff-l). (17) 

The quantity A'off is mainly determined by the nanopore 
length and can be approximated using Eq. (fT5|) . The 
quantity 70 is governed by the specifics of the MC dy- 
namics algorithm. In practice, it can be determined from 
MC measurements of the time-dependence of (Az^) for 
polymers of different lengths confined to an infinite cylin- 
drical tube. 



IV. METHODS 

A. Free Energy Calculations 

Monte Carlo simulations employing the Metropolis 
algorithm and the self-consistent multiple histogram 
(SCMH) method'*'* were used to calculate the free en- 
ergy functions for the polymer-nanopore model described 
in Section nil The results of these calculations were pre- 
sented in a recent article,^ to which the reader should 
refer for a detailed description of the methodology. 

B. MC Dynamics Simulations 

Monte Carlo dynamics simulations were used to study 
the translocation dynamics for the system. Polymer 
motion was generated through random monomer dis- 
placement, in which the coordinates of a randomly cho- 
sen monomer were displaced by an amount Ar> for 
A = x,y, z. Each coordinate displacement was randomly 
chosen from a uniform distribution [— Amax, Amax]- Un- 
less otherwise stated, A^ax = 0.14cr. A trial moves was 
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rejected if it lead to overlap with another nonbonded 
monomer, violation of the bonding constraint, or over- 
lap with the nanopore or barrier wall. The polymer co- 
ordinates were chosen to correspond a desired initial co- 
ordinate, Qq, and the system was equilibrated for fixed 
Qo for typically 10^ MC cycles. Two different first pas- 
sage times were measured: (1) ri, the first time that all 
monomers had completely emptied from either the cis or 
trans domain, and (2) t, the first time that all monomers 
had completely emptied the nanopore and were all in ei- 
ther the cis or trans region. We define T2 to be the time 
taken by the polymer to empty the nanopore, and thus, 
T = Ti + T2. We simulated typically 500 - 2000 translo- 
cation events to calculate mean translocation times and 
lO'' — 10^ events in cases where translocation time distri- 
butions were desired. 



C. Theoretical Calculations 

MC dynamics simulations of a hard-sphere chain poly- 
mer in an infinitely long cylinder were used to measure 
(Az^) as a function of time for the centre of mass of 
a polymer for several lengths. A linear fit to each of 
these functions was used to extract the friction coeffi- 
cient 72 using Eq. (|12|) . Fitting 7^ vs polymer length 
using Eq. ([T^ yielded the constant 70 for a tube of a 
specified radius. 

Equations ([7]), ([5]) and © were used to calculate the 
mean first passage time for the translocation events. The 
integrals were calculated numerically using the trape- 
zoid rule. The friction coefficient 7 was calculated using 
Eq. p7|) . and the corresponding diffusion coefficient, D, 
was then calculated using Eq. p6)) . The value of the ef- 
fective number of monomers in the nanopore, iVeff , was 
chosen so that the calculated (ri) was equal to the value 
measured directly in the MC dynamics simulations. Af- 
ter determining D, Eq. ([5]) was solved numerically using 
a simple finite-difference method with the boundary con- 
dition Q — Qo a.t t = to yield the function 'P{Q,t). 
The distribution of first passage times, gij; Qq) was cal- 
culated by solving Eq. the integral was first solved 
numerically using the trapezoid rule, and the derivative 
was calculated using the finite-difference approximation. 

In the results presented below, lengths are measured 
in units of cr, the monomer diameter. In addition, time 
is measured in MC cycles, where 1 MC cycle corresponds 
to one attempted move per monomer, on average. 

V. RESULTS 

The focus of this study is the elucidation of polymer 
translocation dynamics in the quasi-static regime. In this 
regime, the FP formalism described in the previous sec- 
tion is valid and can be used to calculate first passage 
times associated with translocation. For a given model 
system, it is not clear a priori which parameter values 



correspond to this regime. One effective means to de- 
lineate the regime boundaries is to measure the scaling 
exponent a, defined by the power law scaling relation: 

(T)cxiV". (18) 

In the limit of sufficiently strong pore friction, where the 
system is in the quasi-static regime, it is expected that 
a — 2 for unforced translocation. For a nanopore of finite 
length L, the most appropriate scaling relation is slightly 
different. For a narrow nanopore where the polymer has a 
linear configuration inside the nanopore, the appropriate 
relation is: 

(r) cx (iV - A^p)", (19) 

where 

A^p = L/<T + 1 (20) 

is approximately the number of monomers that lie in- 
side the nanopore. The rationale for this modification is 
that most of the time for translocation passes during the 
phase in which the nanopore is filled, and monomers lie 
outside in both the cis and trans regions, while the time 
elapsed during the pore-emptying stage is very small, by 
comparison. Thus, translocation is almost complete by 
the time either the cis or trans region first empties at 
time Ti, and t = ti + T2 ^ ti, since T2 ^ ti. Justifi- 
cation for this claim will be presented below. From this 
perspective, the effective length of the polymer is N — Np. 
Of course, in the limit that A^ ^ Np, Eq. (fT9)) reduces to 
Eq. (UHl). 

It is instructive to first compare the scaling results in 
this MC dynamics study with those from the Langevin 
dynamics study of Ref. [ttI which examined polymer 
translocation through a hole in a fiat two dimensional 
barrier. Because of the variability of bond lengths per- 
mitted by Eq. we choose a small, but finite, barrier 
width of L = 0.1. This precludes the unphysical pos- 
sibility that two bonded monomers are bisected by the 
wall away from the pore. This barrier width should be 
small enough to approximate a 2-D barrier. Figure [2] 
shows the mean translocation time (r) vs polymer length 
A^ for several different pore radii. In these simulations, 
the polymer was initially placed midway in the nanopore 
such that Qo = 0.5. Over the range of A^ considered 
(A^ = 20-100) (r) satisfies the power relation of Eq. (fTO|) . 
The solid curves show the best fits to the simulation data. 
The inset of the figure shows the corresponding values of 
the scaling exponents, a, as a function of pore radius, 
R. The exponent has a value of a « 2.4 in the range 
R = 0.65-0.9, but decreases rapidly with decreasing R 
for R < 0.65 and reaches a = 2.07 ± 0.01 for R = 0.525. 
These data are in general agreement with the results of 
Ref. [73. Note that this scaling behaviour cannot be ex- 
trapolated to very large A. As A^ increases, the scaling 
is expected to change in a manner such that (t) increases 
monotonically with decreasing R. 

Figure [3] shows similar scaling data for a nanopore of 
length L = 2.917. As in the case described above, de- 
creasing R results in an increase in (t) . Consequently, the 
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FIG. 2. Mean translocation time, (r), vs polymer length, TV, 
for a pore of length L — 0.1 for several different pore radii. 
The inset shows the scaling exponent, a, vs pore radius R 
obtained from fits to the data. The solid red curve is a guide 
for the eye. 



polymer can more easily maintain a state of conforma- 
tional equilibrium as translocation proceeds, which leads 
to the observed decrease in the scaling exponent with R. 
The scaling exponent reaches a value of a = 2.003 ±0.004 
for a nanopore radius of i? = 0.55, which corresponds to 
the case of quasi-static dynamics and strong pore friction. 
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FIG. 3. Mean translocation time, (r), vs N ~ Np, for a pore 
of length L — 2.917 for several different pore radii. The inset 
shows the scaling exponent, a, vs pore radius R obtained from 
fits to the data. The solid red curve is a guide for the eye. 

One reason that decreasing the pore radius increases 
(r) is related to the variation of the friction parameter 
7o with R. As R decreases, the acceptance probability 
of trial displacements of monomers inside the pore de- 
creases because of an increasing likelihood that the move 
leads to overlap with the nanopore wall. This reduces 
the mobility for these monomers and therefore increases 
in 70. This trend is illustrated in Fig. H) From Eq. ([TTl) . 
it follows that the friction coefficient for the translocat- 
ing polymer, 7, should also increase. The inset of Fig. |3] 
shows the variation of the friction coefRcient for a poly- 
mer in an infinite cylindrical tube, 7^, as a function of 
iV — 1 for i? = 0.55, where each 7^ value was deter- 
mined by measuring (Az^) as a function of time and us- 



ing Eq. (IT2I) . Fitting the data using Eq. p3t yields the 
value of 7o. These results confirm the scaling relation of 
Eq. p^ . [Note that the coupling between 70 and i? is a 
consequence of the MC dynamics algorithm used. An al- 
ternative choice of trial moves for monomers in the pore, 
e.g. smaller displacements in the x-y direction compared 
to those along z, would allow decoupling of friction and 
pore radius.] 
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FIG. 4. Friction parameter, 70, vs nanopore radius, R, for 
a polymer in an infinitely long cylindrical tube. The inset 
shows the friction coefficient 7^, vs polymer length, A'^ — 1, for 
a polymer in an infinite cylindrical tube of radius R — 0.55. 
The red curve is a fit to the linear function of Eq. (|17p 

A second factor contributing to slower translocation as 
R decreases can be understood from the effect of decreas- 
ing pore radius on the free energy profiles. As discussed 
in Ref. |4^, the free energy is characterized by an oscil- 
lation with an amplitude, Ai^osc, that increases with de- 
creasing R. This is illustrated in Fig. [5ja), which shows 
profiles for a TV = 20 polymer in a pore of length L = 6 
for pore radii of i? = 0.55 and R = 0.8. The origin 
of the oscillations is related to the nature of the vari- 
ation with Q of the orientational entropy of the bonds 
at the two edges of the nanopore for this model. When 
the oscillation amplitude is of order fee 7" or larger, the 
translocation rate will be appreciably reduced. Thus, de- 
creasing R increases Ai^osc, which in turn increases (r). 
In addition to increasing the translocation times, the free 
energy oscillations can have other effects on the translo- 
cation dynamics. Figure [SJb) shows the history of Q for 
a single translocation event for a A'^ = 20 polymer in a 
pore of length L = 6 and radius R = 0.55. The polymer 
clearly undergoes a diffusive type of motion while it is 
in the plateau region of the free energy profile in Fig. [S] 
Note, however, that the polymer tends to dwell longer at 
Q values close to the local free energy minima and jump 
rapidly over the free energy barriers. As R decreases, 
the oscillation amplitude increases, and the increase in 
the roughness of the landscape decreases the transloca- 
tion rate. 

In addition to its dependence on pore radius, the am- 
plitude of the free energy oscillations also displays an 
interesting dependence on nanopore length. As is evi- 




time 

FIG. 5. (a) Translocation free energy profiles for a polymer 
of length A'^ = 20 in a nanopore of length L — 6 for two 
different pore radii, (b) Q, vs time for an escape event for 
R = 0.55. The horizontal dashed lines mark the values of Q 
corresponding to local free energy minima. 



dent in the inset of Fig. [51 the amplitude itself oscillates 
with L and has a period of AL = 1, i.e. one bond length. 
The origin of this effect was explained in Ref. H^. For 
the case of i? = 0.55 illustrated in the figure, the os- 
cillation amplitude varies from negligible magnitude to 
a maximum of /JAFqsc ~ 2.1. The effect on the varia- 
tion of (ti) and (r) with L is shown in the figure and 
follows the expected pattern: for high values of Ai^osc, 
translocation proceeds more slowly and the translocation 
times are high, while for nanopore lengths corresponding 
to low AFosc, the translocation times are also low. This 
behaviour is expected to be independent of N, as long as 
the pore friction is sufficiently high. Results for N — 100 
with L = 2.917 and L = 2.333 are consistent with this 
pattern (data not shown). 

As explained in Section IIII[ application of the FP for- 
malism requires the value of the translocation diffusion 
coefficient or, equivalently, the friction coefficient, 7, to 
predict translocation times. Since the friction coefficient 
cannot easily be determined, we employ the procedure 
described earlier. The value of 7 is chosen such that the 
mean translocation time, (ti), calculated using Eq. ([7]) 
is equal to the value obtained directly from the dynam- 
ics simulations. This value can then be used in Eq. p7)) 
to calculate N^s, whose value should be close to the es- 
timate in Eq. ([T5|). Deviations from this estimate can 
be interpreted as a breakdown in the quasi-static con- 
dition and/or the assumption that pore friction is the 
dominant contribution to the total friction. The scaling 



2 3 4 5 6 7 

L 

FIG. 6. Mean translocation times, (ri) and (t), vs nanopore 
length, L, for a polymer of length A'^ = 20 and a nanopore 
of radius R = 0.55. The inset shows the variation of the free 
energy oscillation amplitude, /3A_Fosc vs L. 

7 oc A^^ predicted by Eq. ([T7]) provides another such test 
of quasi-equilibrium. Finally, this value of 7 obtained by 
fixing (ri) can be used in Eqs. ([5]) and © to predict the 
complete distribution of times, g{T; Qo). 

Figure [7] shows the variation of the calculated A'cff with 
R for a = 20 polymer in a pore of length L = 3. At 
the high end of the range of values of R shown, A"off is 
significantly larger than the value of A'off ~ 5 predicted 
by Eq. (jl5l) . The average number of monomers inside the 
pore was observed to increase only slightly with R, and so 
it does not account account for this result. Consequently, 
it must be due to decrease in 70 with R, illustrated in 
Fig. m Generally, N^s decreases with decreasing R for 
R > 0.56. Below this value, it is approximately constant, 
with a value A'eff ~ 5.1. This suggests that the quasi- 
static condition holds in this low-i? regime. This inter- 
pretation is consistent with the variation of the scaling 
exponent a with R in Fig. [31 where the quasi-static limit 
value of a = 2 is observed for R = 0.55. The inset shows 
the variation of Aeff with L for the same A^ = 20 polymer 
and a pore of radius R = 0.55. While N^e does tend to 
increase with L, there is a superimposed oscillation with 
a period of AL = 1. To help understand these results, 
the variation of Acff with L was fit to the function 

TVoff = L/o- + ao - ai sin(27rL/cr). (21) 

This function provides a reasonable fit over the range 
of L considered, which yielded values of oq = 2.18 and 
ai = 0.564. Ignoring the oscillation, the prediction of 
A"off = _L/(T + 2.18is roughly consistent with the pre- 
diction of Eq. (IT5|) . The origin of the oscillation is not 
completely understood at present but it may be a result 
of the MC algorithm employed. Note that the values of 
ao and ai are such that A'cff > L/a + 1.62. The lower 
limit exceeds the average number of monomers inside the 
nanopore. Thus, the oscillation in A^eff is associated with 
the reduced mobility of monomers that lie outside the 
pore, the most significant of which lie in the zones of 
partial confinement near the pore edges {—a/2 < z < 
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and L < z < L + a/2). This effect is similar to that 
of the periodic variation of Ai^osc with L illustrated in 
Fig. [6l which is also associated with the monomers near 
the edges 4^ 



7.5 



7 - 
6.5- 

5.5- 
5- 
4.5L- 



7 - 

6 - 



I ' I ' I ' I ' I 




I ■ I ■ I ■ I ■ I 



2.5 3 3.5 4 4.5 5 




0.55 



0.6 



0.65 
R 



0.7 



0.75 



0.8 



FIG. 7. Effective number of monomers in the nanopore, A'cff , 
vs pore radius for a polymer of length A'^ = 20 in a nanopore 
of length L = 3. The inset shows the variation of A'efr with 
L for the same A'^ = 20 polymer in a pore of radius = 0.55. 
The red curve is a fit to Eq. (f2T|) . 
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FIG. 8. Translocation time distributions for a polymer of 
length A'^ = 20 and a nanopore of radius R = 0.55 and length 
L = 3. Distributions for two different Qo values are shown. 
The solid curves are theoretical predictions using A'cff = 5.1, 
which is obtained by setting the mean predicted time to be 
equal to the mean time from the dynamics simulations. The 
inset shows the free energy function for this system. 



consequently, 



Figure [5] shows translocation distributions calculated 
for a iV = 20 polymer and a nanopore of dimensions 
L = 3 and R = 0.55. Distributions for Qo = 0.3 and 
Qo — 0.5 arc shown. The hgure also shows theoreti- 
cal predictions for the distributions, calculated using the 
method described in Sections HIT] and IIVCI The fric- 
tion coefficient 7 was determined using Eq. (1171) with 
A'efi = 5.1, the value obtained by constraining the pre- 
dicted (ti) to be equal to that measured directly from the 
MC dynamics simulations for Qo = 0.5. The theoretical 
distributions are in excellent quantitative agreement with 
those obtained by simulation. 

The translocation time distributions in Fig.|S]each dis- 
play exponential tails at longer times. The origin of 
this feature can be understood in the context of simple 
diffusion between absorbing boundaries. For unforced 
translocation, the polymer spends most of its time mov- 
ing over a relatively flat plateau of the free energy func- 
tion, which is shown in the inset of the figure. When it 
reaches the plateau edge, the polymer quickly exits the 
pore during the pore-emptying stage, i.e, r = T1+T2 ~ n, 
as described above. Neglecting the oscillations and the 
very slight overall negative curvature of the plateau, this 
process can be viewed approximately as simple diffusion 
with absorbing boundary conditions at Qa = 0.1 and 
Qb = 0.9. In this case, Eqs. (O and ^ can be solved 
analytically, yielding^ 



9{t;Qo) = C(Qo)exp(-T/rc), 



(23) 



where 



C(Qo) 



dTrDoff . 
— 72— 



MiQo~Q.)/LQ), (24) 



and where the time constant, Tc, is given by 



(25) 



Thus, the exponential decay time is predicted to be inde- 
pendent of Qo, which is consistent with the data in the 
figure. Fitting the data in the exponential region to the 
function 



(26) 



5(r,Qo) = ^E/^^[l 

^ rj— 1 



where /3p = -kp/Lq, and Lq = Qa — Qh = 0.8. For long 
first-passage time t, the p = 1 term becomes dominant; 



5(T) = C'exp(-r/T^) 

yields « 1.98 x 10'^. Using this value in Eq. (^51) . 
we estimate Z^cff ~ 3.2 x 10^^. Of course, this ap- 
proximation is strictly valid for a perfectly flat free en- 
ergy curve. The corrugated structure of the free energy 
plateau slows down the rate of diffusion, and thus this 
value is expected to be lower than the value of D cal- 
culated using the methods described in Section IIIIl Us- 
ing Eqs. and (IT71) with a value A'cff = 5.1 yields a 
value of D = 4.84 x 10~*, which is greater than Deff, 
as expected. The fit also yields C = 6.7 x 10~^ for 
Qo = 0.5, and C = 4.5 x 10"'^ for Qo = 0.3. These 
values are approximately consistent with the values pre- 
dicted using Eq. ^ of C(Qo = 0.5) = 6.3 x 10"^ and 
HSpLq)] sin(/3p(Qo - Qa)) expf^^m^^} = 4.4 x 10"^. The smafl discrepancies arise 

from the approximations employed. 

For the distributions shown in Fig.[Sl it is evident that 
the mean translocation time for Qo = 0.3 is lower than 
the mean time for Qo = 0.5. This is an illustration of a 
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general result: the closer the polymer is initially situated 
to either one of the pore exits, the shorter the translo- 
cation time. This arises because the polymer exits more 
frequently at the closer pore boundary, and the shorter 
distance reduces the first passage time. The longer time 
required to reach the more distant exit is offset by the 
diminished probability of reaching that exit. This is con- 
sistent with the predictions for simple diffusion between 
absorbing boundariesji It also effects a shift of the entire 
distribution to shorter times. 

Figure. [H]shows the variation with Qq of the probability 
of a trans side exit. As the initial position is moved closer 
to the cis side, i.e. as Qq decreases, the probability of a 
trans exit decreases, and the cis exit probability increases 
accordingly. The theoretical curve was calculated using 
Eq. with Qa = 0.1 and Qh = 0.9, the bounds of 

the free energy plateau shown in the inset of Fig. [51 The 
theoretical prediction is in excellent agreement with the 
simulation probabilities. The observed variation of TTtrans 
with initial position corresponds closely to the prediction 
for a diffusive process with a constant free energy between 
two absorbing boundaries, where the variation is exactly 
linearji The slight deviation from linearity observed in 
Fig. [S] arises from the overall slight negative curvature of 
the (rough) free energy landscape evident in the inset in 
Fig.lH 




FIG. 9. Translocation probability, TTtrans vs initial coordinate 
Qo for a polymer of length N = 20 and a nanopore of radius 
R = 0.55 and length L = 3. The symbols are data calculated 
directly from the MC dynamics simulations, and the solid line 
is a theoretical prediction using Eq. H10|) . 

Figure [TO] shows translocation time distributions calcu- 
lated for three different polymer lengths using a nanopore 
with L = 3 and R = 0.55 and an initial coordinate of 
Qo — 0.5. As the polymer length increases, the distri- 
bution tends to broaden and the mean value increases. 
The theoretical distributions overlaid on the simulation 
data were all calculated using a friction coefficient, 7, 
obtained using Eq. (|I7I) with N^ff = 5.1. As explained 
in Section IIIIl the value of A^eff is expected to depend on 
the nanopore dimensions, but not on the polymer length. 
The excellent agreement between the theoretical distribu- 



tions and those calculated by MC dynamics simulations 
using the same N^ff is consistent with this expectation. 




5 10 15 20 25 

T/IO' 



FIG. 10. Translocation time distributions for polymers of 
lengths A'^ — 20, 26 and 30, and a nanopore of radius _R = 0.55 
and length L — 3. The solid curves are theoretical predictions 
using A^efi = 5.1. The inset shows the free energy functions for 
the three systems. The functions have been vertically shifted 
to prevent overlap for visual clarity. 

Figure [TT] shows a comparison of translocation distri- 
bution times for nanopore lengths of L = 3 and L = 2.5, 
for the case where the polymer length is A^ = 20 and the 
nanopore radius is R = 0.55. The distribution is shifted 
to lower times for the shorter nanopore. The theoretical 
curves shown in the figure show excellent agreement with 
the simulation data. Note that the theoretical curve for 
L — 2.5 was calculated using Aoff = 4.6, a value that is 
consistent with the estimate in Eq. . A fit of the ex- 
ponential tail of the L = 2.5 distribution yielded a time 
constant of = 1.135 x 10^. As seen in the figure in- 
set, the free energy function is much smoother over the 
plateau region for L — 2.5 than for i = 3. In addi- 
tion, the width of the plateau is approximately the same, 
i.e. Lq « 0.8 in each case. Employing the model of 
free diffusion between absorbing boundaries at Q — 0.1 
and 0.9 for the case of L = 2.5, Eq. ([^5)1 yields a value 
of Dcs = 5.7 X 10~*, which is very close to the value 
D ^ 5.48 X 10"*^ predicted using Eqs. I|16[) and (|17p . 
This is exactly the consistency expected for a smooth, 
flat free energy function, and stands in contrast to the 
discrepancy between similarly calculated values of D^s 
and D for L = 3 discussed above in the context of Fig. [5] 
It follows that much of the reduction in the translocation 
time upon decreasing the nanopore length from L = S to 
L = 2.5 is accounted for by decreasing roughness of the 
free energy function; a smaller contribution arises from 
the small decrease in A'cff and corresponding small in- 
crease in D. 

Figure [T^j shows translocation time distributions for 
a polymer that translocates out of a spherical cavity, 
the system illustrated in Fig. [Ijb). The cavity radius 
is Rc = 3, and as before, A^ = 20, L = 3, and R = 0.55. 
Distributions for three different Qq are shown. Compar- 
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FIG. 11. Translocation time distributions for a nanopore of 
radius R — 0.55 and length L = 3. The polymer started 
midway through the pore with Qo = 0.5. Distributions for 
two different nanopore lengths are shown. The solid curves 
are theoretical predictions using A'^ofi = 5.1 for L = 3.0 and 
iVcff = 4.6 for L = 2.5. 



ing the results with those of Fig. [5] for the flat barrier 
system shows that translocation times are slightly lower 
in the spherical cavity system. The origin of this effect 
is clear from a comparision of the free energy functions, 
shown in the figure inset. The confinement of the poly- 
mer in the cavity breaks the symmetry and leads to a 
free energy gradient. The resulting entropic driving force 
pushes the polymer to the trans side and reduces the time 
taken for the polymer to exit the nanopore. 




10 



x/10" 



FIG. 12. Translocation time distributions for the spherical 
cavity system as illustrated in Fig. [T^b) for a cavity radius of 
Rc = 3. In addition, R = 0.55, L = 3 and iV = 20. The theo- 
retical curves were calculated as described in the text, using 
NcS = 5.6. The inset shows the free energy function used in 
these theoretical calculations as well that used in calculations 
for the planar geometry. 

The asymmetry in the free energy function means that 
Qo — 0.3 and Qo = 0.7 no longer have the same distribu- 
tions, unlike the case for the symmetric flat-barrier sys- 
tem of Fig. [Tfa). Instead, there is a monotonic decrease 
in the distribution mean times as the starting point of the 



polymer is placed closer to the trans side of the nanopore 
(i.e. larger Qo)- The bias introduced by the entropic 
driving force causes the majority of translocation events 
associated with a polymer exit on the trans side of the 
pore. The closer to the trans side the polymer starts, 
the faster will be its exit. Note that this trend must in- 
evitably reverse for sufficiently low Qo: i-e. if the polymer 
is almost complete inside the cis cavity, in which case the 
polymer may exit very rapidly into the spherical cavity. 

The theoretical distributions shown in Fig. [T^] are in 
excellent agreement with the simulation data. Note that 
the theoretical distribution for the spherical cavity sys- 
tem shown in the figure was calculated using A'eff = 5.6, 
somewhat larger than the value of N^ff = 5.1 employed 
for the flat-barrier system. As discussed above, iVeff 
is expected to depend only on L, which has the same 
value {L — 3) as that used in the flat barrier simula- 
tions. The slightly higher value of iVoff could be a result 
of out-of-equilibrium effects which arise from more rapid 
translocation due to the entropic driving force, as well as 
an increase in the relaxation time of the section of the 
polymer inside the spherical cavity. If so, then D and 
7 would likely exhibit a dependence on polymer length. 
As a check, we calculated the translocation time distri- 
butions for a longer polymer, N — 30, for an otherwise 
identical system. The theoretical predictions are in com- 
parably good quantitative agreement with the distribu- 
tions calculated from the dynamics simulations for the 
same value of Ncff (data not shown). This suggests that 
the difference in the values of TVeff for the two systems is 
not an out-of-equilibrium effect. More likely, the effect 
arises from the curvature of the cavity wall in the vicin- 
ity of the cis entrance to the nanopore. The curvature 
should lead to a slightly more crowded region near the 
entrance, which will lead to a slightly lower acceptance 
rate in the MC algorithm due to an increased likelihood 
of particle overlap. This in turn leads to a decrease in the 
monomer mobility near the entrance and consequently a 
slightly larger 7 and Ngg. It has already been noted in 
the discussion above regarding the data in Fig. [7] how 
sensitive the pore friction can be to very small changes 
in the nanopore characteristics. 

The FP predictions of translocation time distributions 
are in excellent quantitative agreement with those mea- 
sured using MC dynamics method for the hard-potential 
model described in Section |TT1 In order to demonstrate 
the generality of this result, we present preliminary re- 
sults for a repulsive Lennard- Jones chain model. The de- 
tails of the model are essentially the same as in Ref. [ttI . 
except that nanopore of finite length is used. The free 
energy function used in the FP calculation was calcu- 
lated as in Ref. '46^, and Brownian dynamics simulations 
were used to measure the distributions. Quasi-static con- 
ditions were achieved by increasing the friction coeffi- 
cient of monomers inside the pore, 7p, relative to that 
for monomers ouside, 70. Figure [T3] shows results for 
7p/7o = 16, for which there is excellent quantitative 
agreement between theory and simulation. The inset 
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shows that the measured Ncs approaches the expected 
value for sufficiently large jp/jo- The results demon- 
strate that the FP methodology provides a valid descrip- 
tion of quasi-static translocation dynamics for this com- 
monly used simulation model. Further details will be 
presented in a future study. 
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FIG. 13. Translocation time distribution for a repulsive LJ 
chain for the nanopore geometry of Fig. [Tfa) calculated using 
Brownian dynamics simulations. Here, 7p/7o = 16, i? = 0.65, 
i = 3, and A'' = 20. The solid line is the theoretical curve 
calculated using the FP formalism. The inset shows the vari- 
ation of A'efi with 7p/7o. 



VI. DISCUSSION 

In this work, we have studied the translocation dy- 
namics of a simple polymer/nanopore model system us- 
ing Monte Carlo dynamics simulations and the Fokker- 
Planck formalism. The emphasis was on the quasi-static 
regime, in which the FP formalism provides a valid ap- 
proach for describing translocation. While many theo- 
retical studies have employed this approach, almost all 
have used approximate analytical forms for the free en- 
ergy function. In our study, we use free energy func- 
tions calculated directly using a MC multiple histogram 
method for the same model used in the dynamics simu- 
lations. The translocation time distributions calculated 
using these functions thus provide a much better test of 
this theoretical approach than would otherwise be possi- 
ble. 

For the MC dynamics algorithm employed here, the 
key parameter governing the dynamical regime is the 
nanopore radius, R. As noted in a previous study,^- de- 
creasing R can slow the translocation rate in a manner 
that leads to a decrease in the translocation time scaling 
exponent, a. In the present case, this happens in part by 
effectively increasing the friction of the polymer in the 
pore. For sufficiently low R, the translocation rate de- 
creased to the point that a value of a = 2 was observed, 
indicating that the quasi-static regime was attained. This 
effect is similar to that observed in Ref. ISll in which 



decreasing the solvent viscosity was used to attain the 
same exponent value. In each case, quasi-equilibrium is 
observed when the ratio of the translocation time to the 
conformational relaxation times of the cis and trans chain 
sections is sufficiently high. Muthukumar has noted that 
this is the relevant dynamical regime for experiments of 
DNA translocation through a-hemolysin pores. ^ 

Knowledge of the translocation free energy functions 
can provide valuable insight into the quantitative trends 
of the translocation dynamics. An useful illustration in 
this study is the effect of oscillations in the functions 
that occur for sufficiently small pore radius, The 
translocation time is generally longer for systems with 
a higher free energy oscillation amplitude. Interestingly, 
the oscillation amplitude itself can be very sensitive to 
the nanopore length, L, as is evident in Fig. [51 which in 
turn leads to a sensitive dependence of the translocation 
time on L. While it is unlikely that oscillations arising 
from the same physical origin are present in the systems 
studied experimentally^ this effect may be highly rel- 
evant to interpretation of results from other simulation 
studies, many of which use qualitatively similar models 
to that employed here. Since the oscillations are asso- 
ciated with the orientational freedom of only two poly- 
mer segments localized at the nanopore edges, the effect 
of slowing the translocation rate may even be present 
in cases where the trans and/or cis chain sections are 
out of equilibrium. Consequently, small changes in the 
nanopore geometry (such as L, here) could lead to sig- 
nificant changes in the translocation times, with corre- 
sponding changes in the scaling exponent, even outside 
the quasi-static regime. Such effects may account in part 
for the wide range of scaling exponents reported in vari- 
ous simulation studies. This underscores the importance 
of calculating accurate free energy functions to aid in the 
interpretation of dynamics simulation data. 

Translocation time distributions for systems deter- 
mined to be in the quasi-static regime were calculated 
using the FP formalism with measured free energy func- 
tions. The calculations assume that the nanopore friction 
is the dominant contribution to the overall friction coef- 
ficient for the translocation process. The distributions 
were calculated using a coefficient 7 which was charac- 
terized by the parameter iVoff, the effective number of 
monomers whose motion is strongly affected by nanopore 
confinement. In the case of unforced translocation, we 
investigated the effects of varying polymer length, N, 
nanopore length, L, and the initial position of the poly- 
mer, Qq. In addition, we measured distributions for poly- 
mers translocating out of a spherically confined space, a 
feature that gives rise to an entropic force that drives 
translocation. The calculated distributions were gener- 
ally in excellent quantitative agreement with those ob- 
tained directly from dynamics simulations for physically 
meaningful values of iVeff- Preliminary results of calcula- 
tions using a repulsive LJ chain chain model and Brown- 
ian dynamics simulations yielded comparably good agree- 
ment. 
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In this study, calculations were carried out mainly for 
short polymer chains. This was due to the consider- 
able computational effort required to obtain very accu- 
rate free energy functions and translocation time distri- 
butions from the simulations. In future work, we will in- 
vestigate the effectiveness of FP formalism in describing 
translocation dynamics of much longer polymers. An- 
other important issue not adequately addressed in this 
study is the effect of a driving force, which plays an es- 
sential role in DNA translocation experiments and whose 
effects have been extensively investigated in other simu- 
lation studies. Though we have briefly considered the 
effect of an entropic driving force in this work, a driving 
force associated with an potential energy gradient across 
the nanopore is not suitable for the MC dynamics algo- 
rithm employed here.^^ In a future study, we will examine 
this effect using Langevin and Brownian dynamics sim- 
ulations. In addition, we will more thoroughly elucidate 
the regime within which the quasi-static approximation is 
sufficiently valid to justify using the FP formalism. This 
regime is expected to show a complex dependence on the 
polymer length, the conformational relaxation times of 
the chain, pore friction, the driving force strength, as 
well subtle features of the free energy functions. This 
work will also be extended to using more realistic molec- 
ular models to provide a more meaningful comparison 
with experimental systems. The long term goal of this 
research is to provide accurate predictions of transloca- 
tion times using precisely calculated free energy functions 
within a well defined regime of validity. 
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